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We investigate the statistics of extremal path(s) (both the shortest and the longest) from the root 
to the bottom of a Cayley tree. The lengths of the edges are assumed to be independent identically 
distributed random variables drawn from a distribution p(l). Besides, the number of branches from 
any node is also random. Exact results are derived for arbitrary distribution p(l). In particular, for 
the binary {0, 1} distribution p(l) = p<5i,i + (1 — p)Si,o, we show that as p increases, the minimal 
length undergoes an unbinding transition from a 'localized' phase to a 'moving' phase at the critical 
value, p = p c = 1 — 6 _1 , where b is the average branch number of the tree. As the height n of the 
tree increases, the minimal length saturates to a finite constant in the localized phase (p < p c ), but 
increases linearly as v m i n (p)n in the moving phase (p > p c ) where the velocity i; m j n (p) is determined 
via a front selection mechanism. At p = p c , the minimal length grows with n in an extremely 
slow double logarithmic fashion. The length of the maximal path, on the other hand, increases 
linearly as i> ma x(p)n for all p. The maximal and minimal velocities satisfy a general duality relation, 
?Wn(p) + «max(l — p) = 1, which is also valid for directed paths on finite-dimensional lattices. 



PACS numbers: 05.40.-a, 64.60. Cn, 02.50.-r 
I. INTRODUCTION 

Optimization lies at the heart of a vast number of 
phenomena: Particles "seek" paths with minimal action, 
species "try" to maximize fitness, companies minimize 
cost. Many optimization problems, e.g., the traveling 
salesman problem, are notoriously hard . Additionally, 
optimization problems often involve randomness which 
makes them even more complicated. When the number 
of random entries is very large, however, optimization 
problems might become simpler thanks to the "concen- 
tration of the measure" . This means that the probability 
distribution of some random variable becomes highly lo- 
calized, almost like a delta function concentrated around 
its average value. This phenomenon is well-known in 
probability theory where, for instance, it accounts for 
the fact that after flipping a coin N times, the probabil- 
ity that the number of heads exceeds N/2 by more than, 
say, lOOV^/V is about 10~ 800 . Below, we shall investigate 
an optimization problem in random media which demon- 
strates this concentration of the measure in the strongest 
sense (the variance of the optimized quantity is finite) 
and displays an unbinding phase transition. A remark- 
able hidden connection with traveling wave phenomena 
allows to explain both the concentration of the measure 
and the phase transition. 

The optimization problem considered in this work can 
be formulated as follows: Take a rooted tree whose bonds 
have random lengths and find descending paths of ex- 
tremal total length. We assume that a random length 
(or energy) is assigned to each bond of the tree, and that 
lengths are independent and chosen from the same prob- 
ability distribution p(l). The total length of a path is the 
sum of lengths of the bonds along the path, and we want 



to determine the minimal (maximal) length among the 
paths from the root to the bottom of the tree. We focus 
on a rooted Cayley tree where each node, except the root, 
has coordination number 3; the coordination number of 
the root is 2. However, our main results can be easily 
generalized to the case of a tree with arbitrary coordina- 
tion number, and even to trees where the coordination 
number at any node is random. 

Similar problems arise in numerous fields ranging from 
computer science Q to condensed matter physics 
where it is known as the problem of directed polymers 
in random media. Indeed, the minimal path problem 
can be considered as the zero temperature limit of the 
problem of directed polymer on Cayley trees. Directed 
polymers on Cayley trees have been investigated in the 
past m, and recently this problem has resurfaced in a 
surprisingly large number of apparently unrelated prob- 
lems 0-^). However, the emphasis of Refs. was on 
the spin glass like transition occurring at a finite temper- 
ature. Here, we consider exactly zero temperature and 
focus on the unbinding transition driven by the parame- 
ter p of the bimodal distribution, p(l) = p5i t \ + (1 — p)<5/.o- 
Recently, similar unbinding (or "depinning" ) phase tran- 
sitions have been found in a number of non-equilibrium 
processes without quenched disorder [|lT|~|l4[ . 

Let us first outline our main results. We derive exactly 
the statistics of both the minimal and maximal path on 
a random Cayley tree for arbitrary distribution p{l) of 
the edge lengths. We find a class of distributions for 
which the minimal path undergoes an unbinding phase 
transition from a localized phase to a moving phase as 
a parameter of the distribution is varied. A particular 
example, which we study in detail in this paper, is the 
bimodal distribution. In the bimodal case, length of each 
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bond can be either 1 with probability p or with prob- 
ability (1 — p). For any path starting at the root and 
moving down, the length of the path is the sum of the 
lengths of the bonds along the path, i.e., the number of 
l's along the path. There are 2™ paths from the root 
to the n th level. We show that the minimal length ex- 
hibits an unbinding transition as the parameter p varies 
between and 1. There exists a critical value p c = 1/2 
such that the average minimal length grows linearly with 
n when p > 1/2, but gets pinned, i.e., remains finite for 
large n when p < 1/2. More specifically, we find that for 
large n, the average minimal length (L™ m ) behaves as 



(^r n ) 



(ln2)- 1 rnhin 
finite 



p>l/2, 

P=l/2, 
p< 1/2. 



(1) 



For p > 1/2, we shall later compute the minimal "veloc- 
ity" Vmin(p) exactly. The average maximal length, on the 
other hand, behaves as 



i(p)n 



for all p, 



(2) 



where v max (p) + Wmin(l — p) = 1- This duality relation 
is very general and valid for directed paths on arbitrary 
graphs, i.e., for finite-dimensional lattices. We also de- 
rive a generalized duality relation between minimal and 
maximal velocities for a certain class of bounded distri- 
butions. 

The minimal length problem on a rooted Cayley tree, 
where the coordination number at each node is random, 
also displays an unbinding transition at p c = 1 — b~ , 
where b is the average number of branches per node. At 
the critical point, the average minimal length (£™ m ) ex- 
hibits the same double logarithmic growth for any b > 1, 
and above p c the minimal length grows linearly with n 
with velocity v m - ln (p, b). 

The rest of this paper is organized as follows. In Sec- 
tion II, we study the unbinding transition of the minimal 
path for the special bimodal distribution in detail. In Sec- 
tion III, we study the maximal path for the bimodal case 
and derive the duality relation. In Section IV, we gener- 
alize the results to a Cayley tree with random branching. 
In Section IV, we go beyond the bimodal distribution 
and obtain generalized results for arbitrary distribution 
of the edge length. Finally we conclude in Section VI 
with a brief summary and outlook. 



II. MINIMAL PATH(S) 

Instead of studying P™ hl {x) = Prob(L™ in = x), it 
proves convenient to consider the cumulative distribu- 
tion, P n (x) = Prob(L™ in > x). Clearly, P n (x) is the 
probability that all possible paths to the n th level have 
lengths > x. Once we get P n (x), the distribution of 



the minimal length is found from the obvious identity 
P™ ta (x)=P n (x)-P n (x + l). 

It is easy to see that P n (x) satisfies the following re- 
cursion relation 
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(X) = [(1 - P )P n (x) + P Pn(x - 1)Y 



(3) 



This relation can be derived by analyzing various pos- 
sibilities for the lengths of the two edges issuing from 
the root and taking into account that two subsequent 
daughter trees are statistically independent. For in- 
stance, both lengths are equal to zero with probability 
(1 — p) 2 , and paths in the subsequent trees have lengths 
> x with probability P 2 (x); this provides the contribu- 
tion (1 — p) 2 P 2 {x) to P n+ i(x). Similarly, one gets the 
contributions 2p(l —p)P n (x — l)P n {x) and p 2 P 2 (x — 1), 
which sum up to Eq. (3). 

We must solve Eq. (3) subject to the initial condition 



P 



x < 0, 
x > 0. 



(4) 



Clearly, the minimal length L™ ln is a random variable 
which takes values between and n, and Eqs. (d)-(H) 
indeed show that P n (x) — 1 for x < and P n (x) = for 
x > n. 

When n grows, £™ m increases as well and thus it 
should reach a limit. For p < 1/2, this limit is finite, 
i.e., one can find an infinite path with only finite number 
of edges of length 1. Mathematically, it means that P n (x) 
approaches a stationary distribution, P n {x) — ► P(x), 
which satisfies 



P{x) = [(1 - p)P{x) + P P(x - I)} 2 



(5) 



Starting from P(0) — 1, one computes P(x) recursively: 
= p 2 (l — p)~ 2 , etc. P{x) vanishes extremely fast in 
the large x limit: lnP(x) ~ 2 x lnp. 

For p > 1/2, the minimal length grows as n increases. 
To understand why this is so, recall that P n (0) = 1 and 
look at P„(l). By inserting F„(l) = 1 — q„ into Eq. (||) 
we find q n+ i = 2(1 -p)q n - (1 -p) 2 ql- Hence P„(l) -> 1 
for p > 1/2. Proceeding this line of reasoning one can 
verify that P n {x) — > 1 when n — » oo. Additionally, we see 
the difference between the cases of p > 1/2 and p = 1/2: 
In the former situation, P n (x) — > 1 exponentially fast, 
while in the latter situation the approach is algebraic. 
Therefore, the size of the region where P n (x) is close to 1 
should grow linearly with n when p > 1/2. Thus the dis- 
tribution P n (x) approaches the traveling wave form (see 
Fig. (1)), 



P n (x)^V 2 (y), y = x-(L™»). 



(6) 



In Eq. (||), we have chosen y = x — (£™ m ) as the wave 
variable. One can justify this choice by noting that the 
relation P™ in (x) = P n (x) - P„(x + 1) leads to identity 



X [P n {x) - P n (x + 1)] = J2 P n(3 



(7) 
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and hence {L™ ln ) is indeed an appropriate characteriza- 
tion of the location of the front. By inserting the wave 
form (||) into Eq. (||) and using (L™ m ) ~ vn we find, 



T(y -«) = (!- p)^ 2 ^) + P T 2 (y - 1), 



(8) 



The solution of this equation looks like a [1 — 0] waveform 
with P(y) — > 1 as y — > — oo and P(y) — > as y — > oo. 
While it is very hard to solve this nonlinear, nonlocal 
equation exactly, the velocity can be found by analyzing 
the tail region, y — > — oo, where 1 — V(y) is small. Lin- 
earizing Eq. (^) in this region and noting that it admits 



an exponential solution, 1 — V(y) 



, one finds that 



the velocity v\ is related to the decay exponent A via 



vx 



In [2{l-p) + 2pe- x ] 



A 



(9) 



When A > In 2 2 ^ 1 , the velocity v\ is positive. While 
any such A is in principle allowed, and thus the entire 
velocity range of is feasible, a particular value is actually 
selected. This is similar to the velocity selection in a large 
class of problems with a traveling wave front ||l5|]l6l . It is 
well known that for a wide class of initial conditions, the 
extremum value is chosen. From this general front selec- 
tion principle one can infer that in the present case the 
selected value A — A* (p) corresponds to the maximum of 
v\, and hence the selected "velocity" IS V m [ n (p) = v x *. 
Thus, A* is a solution of the equation 



In 



2(1 -p) + 2pe~ 



pX*e~ 



1 — p + p e A * 



= 0, 



and the selected velocity is 



V m in(p) = 



pe 



1 — p + pe 



(10) 



(11) 



Although it seems impossible to provide an explicit ex- 
pression for v m i n (p), one can easily probe its behavior in 
the limiting cases P [ \ and p f 1. In these respective 
limits, one gets 



2p- 1 
In 



and 



v m in(p) 



2p-l 



In 2 

l-p 



(12) 



(13) 



One can also derive the leading correction to the dom- 
inant linear behavior of (L™ m ), 



i{p)n+— Inn. 



(14) 



The logarithmic correction becomes especially important 
when p J. I , as can be seen from the asymptotics 
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ln(l/e) 2 1n(l/e)' 



2p-l. 



(15) 



The logarithmic correction to the front position was 
first derived by Bramson in the context of a reaction- 
diffusion equation |l7],[l5|], and was subsequently re- 
derived and generalized by a number of authors [|l8] g0| . 
Our derivation of Eq. ([lj]) follows an approach of Ref. 
fill , and we only outline the main steps. First, we need 
to examine the exact equation (||) rather than the n — > oo 
limit, Eq. (||), which was sufficient to determine the ve- 
locity. However, the analysis is still performed in the 
region far behind the front, so we can use the linearized 
version of Eq. (||). Writing P n (x) = 1 — Q n (x), plugging 
it into Eq. (0), and ignoring quadratic terms gives 



Q n+1 (x) = 2(1 -p)Q n (x) + 2pQ n {x - 1). 



Now we seek a solution which has the form 



Q n (x) = n a G(yn- a ) 



(16) 



(17) 



with y = x — (L™ m ) as previously and yet undetermined 
exponent a and the scaling function G{z). The correc- 
tion to the leading term in the front position should be 
incorporated as, 
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Vminip) n + c(n) 



(18) 



where the function c(n) is yet undetermined. Plugging 
Eqs. into Eq. @ we find that different lead- 

ing orders are compatible provided that a — 1/2 and 
c(n) = [3\nn with some constant (3. Additionally, the 
scaling function G(z) satisfies a parabolic cylinder equa- 
tion 

^ + z ^ + (2/3A* - l)G(z) = 0. (19) 

Equation ( |l9| ) should be solved subject to the boundary 
conditions G(z) — > as z — > oo as Q n (x) must vanish for 
x — > oo and G(z) ~ z for z —> to ensure that Q n (x) 
is independent on n in the limit n — > oo. The boundary 
condition in the large z limit selects one of the two possi- 
ble solutions: G(z) = Aexp(— z 2 /4) D 2 (px*-i){ z )i where 
D v is the parabolic cylinder function with index v. The 
second boundary condition G(z) ~ z fixes the index of 
the parabolic cylinder function, 2(/3A* — 1) = 1. Hence, 
(3 = 3/2 A* thus completing the derivation of Eq. (fil|). 
Note that the sign of the logarithmic correction term in 
Eq. (|l4|) is positive. This is different from the most ear- 
lier studied problems ]l^-|20|| although the positive sign 
has been occasionally seen, see Ref. [^lj . 

In the the critical case of p = 1/2, the convergence 
of the distribution P n (x) towards the asymptotic value 
is only algebraic, and therefore the front propagates ex- 
tremely slowly. The simplest way to determine the rate 
of propagation is again to look at P n {x) far behind the 
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front. Writing P n (x) 
Eq. © yields 

Qn+l(x) = Q r 



= 1 — Qn( x ) an d plugging it into 
(x)+Q n (x-l) 



i(x) + Q„(x - 1) 



(20) 



Qn{x) = for x < 0, so the first non-trivial <5„'s are 

Q n (l)=q{n). FromEq. ©, g(n+ 1) = ?(n) - ^2. In 
the large n limit, we can employ the continuum approxi- 
mation to get 4* = — ^- and thence Q ra (l) — 4n~ 1 . Sim- 
ilarly, we find the asymptotics Q n (2) ~ 4n~ 1 / 2 . Gener- 
ally, Eq. (|20| ) implies Q n (x) ~ 2^/ Q„(x — 1), from which 



Q n (x) ~ 4n 



. 2 -(*-i) 



= 4 cxp (— 2e 



In In n— a~ In 2\ 



(21) 



This demonstrates that at the critical point, Q n (x) also 
has a traveling wave form, Q n (x) = g(x — x n ), with the 
front at x n ~ (ln2) _1 lnlnn. 

Equation (|2lj) formally applies for x <C x„. To inves- 
tigate Q„(a;) in the entire range of x, we again use the 
fact that the distribution Q n (x) should approach a trav- 
eling wave form, Q n (x) — > Q(y) with y = n~ 2 ( , as 
n — > oo. By inserting this into (^o|) and taking n — > oo 
limit, we arrive at 



Q(y) = 2v/Q(7) 
From this equation, one finds 



Q(y 2 )- 



Q(y) =4y-4y 2 -2y 3 + ... for y I 0. 



(22) 



(23) 



The first term in this expansion indeed agrees with 
Eq. (^l|). Similarly, one gets 



1 - Q(y) ~ ex P 



const 
T^y~ 



for y | 1. 



(24) 



To provide a more rigorous derivation of the growth 
law for (£,™ in ), note that (L™ in ) = J2 x >i p n{x) = n- 
J2x>i Qn{x). Inserting Q n (x) = Q(y) and replacing the 
sum by an integral (valid for large n), we get, 



In n 



Q(e- Z )-. 

2"- 1 lnn z 



(25) 



In the limit of large n, the most important contribution in 
the integral comes from its lower limit. Using Q(y) — » 1 
as y — > 1 (see Eq. (J2J)), we arrive at the desired result, 
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(ln2) _1 lnlnn. 



(26) 



III. MAXIMAL PATH(S) 

We now turn to the statistics of i™ ax , the length of 
the maximal path for the bimodal distribution. It is now 



convenient to define R n (x) — Prob(L" lax < x). By pro- 
ceeding as in the minimal length problem, it is easy to 
show that R n (x) evolves according to the same equation 
(0) with P replaced by R, i.e., 



Rn+l(x) - [(l-p)R n (x)+pR n (x - l)Y 



(27) 



The only difference from the minimal case is in the initial 
condition. Instead of Eq. (Q), we now have 



x>0, 
x < 0. 



(28) 



This difference in the initial condition, however, leads to 
a different behavior of the maximal front as demonstrated 
below. 

Not surprisingly, equation ( |27| ) admits a traveling wave 
solution, R n {x) = R(y = x — v max n) for all p (see Fig. 
(2)). However, in contrast to the [1 — 0] waveform of the 
minimal case, the function R(y) now looks like a [0 — 1] 
waveform with R(y) — > as a double exponential behind 
the front (y — * — oo) and R{y) — > 1 exponentially ahead 
the front (y — > oo). We then analyze the wavefront near 
the 'forward' tail, as opposed to the 'backward' tail of 
the minimal problem. Substituting 1 — R(y) ~ e~^ y , we 
find that there exists a family of traveling wave solutions 
with velocity t> M parametrized by //, 



ln[2(l-p) + 2pe"] 



(29) 



This velocity i> M has a unique minimum at u = fi*. By 
the front selection mechanism, this minimum value v^* 
is selected as the velocity i> max (p) of the front. Thus, we 
have 



In 



2(1 -p) + 2pe h 



p fj, &■ 



1 — p + pe IJj 
and the selected maximal velocity is 



pe H 



1 — p + pe^" 



0, (30) 



(31) 



Using these results together with Eqs. ( |To| ) and (^Tj), one 
immediately finds 



(l-p) = l. 



(32) 



This is the duality relation. Unlike the behavior of the 
minimal length, the maximal length does not undergo an 
unbinding transition in the sense that it increases linearly 
with n for large n for all < p < 1. However, there is still 
a transition at p — 1/2: The velocity w max (p) increases 
from to 1 as p increases from to 1/2 and then sticks 
to 1 for p > 1/2 (see Fig. (3)). 

For the maximal front also, one can easily compute the 
sub-leading logarithmic correction to the front position. 
Proceeding as in the minimal case, we find 



4 



(33) 



Note the negative sign in the correction term, as opposed 
to the positive sign for the minimal case in Eq. (fTi]) . This 



is also consistent with the duality relation in Eq. (32). 
Thus, while the coefficient 3/2 of the logarithmic correc- 
tion term seems to be universal not just in the present 
problem but in many other cases [|l8|-pcfl , the sign of the 
correction term is not and can either be positive |pl| or 
negative 



The duality relation ( |32[) can also be derived by a gen- 
eral argument which does not depend on the structure 
of the underlying tree and therefore is valid for extremal 
directed paths on arbitrary graphs, in particular, on arbi- 
trary lattices. The argument follows from the observation 
that for bimodal distribution, if one replaces the O's by 
l's and the l's by O's on the minimal path (and thereby 
replacing p by 1 — p) , then the minimal path becomes 
the maximal path. Let uq and ri\ denote respectively 
the number of O's and l's on the minimal path. Evi- 
dently, no + n\ — n. Then clearly, L™ m (p) — n\ and by 
duality, I/™ ax (l — p) — n . Adding the two quantities, we 
get L™ m (p) + L" lax (l — p) = n. Dividing by n, we imme- 
diately get the duality relation Eq. ([32]) . Below, we shall 
derive a more generalized version of this duality relation 
which goes beyond the bimodal distribution. 



IV. GENERALIZATION TO TREES WITH 
RANDOM BRANCHING 

The generalization of the above analysis to the case of 
the rooted Cayley tree with a random number of branches 
is straightforward. Let b m is the probability that the 
number of branches is equal to m and b is the average 
number of branches, b — X)m>i TO ^ m - Equation (^) is 
replaced by 

oo 

P n+1 {x) = & ™ [(! ~P) P n(x) +pPn{x ~ I)]™ , (34) 



and the subsequent analysis repeats the steps detailed 
in the case of the binary tree. In particular, P n (x) ap- 
proaches the stationary distribution for p < p c , while 
above p c the minimal length grows according to Eq. ( |l4] ) 
with A* a solution to equation 
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6(1 — p) + bpe 



-x* 



pX*e 



-X" 



1 — p + p e 



-A* 



(35) 



and v m i a (p) given by the same relation (|ll|). At the crit- 
ical point p c = 1 — the minimal length is given by 
the same expression ( p6| ) for any tree. Also, the velocity 
of the maximal path can be determined from the general 
duality relation in Eq. (|32|). 



V. GENERALIZATION TO ARBITRARY 
DISTRIBUTIONS 

We now consider the statistics of extreme path(s) on 
a binary tree for arbitrary distribution p(l) of the edge 
lengths. Let us first consider the minimal path(s). The 
velocity of the minimal path for an arbitrary distribu- 
tion can be extracted from earlier results by Derrida and 
Spohn H which we re-derive below for the sake of com- 
pleteness. We consider again P n (x) — Prob(L™ m > x). 
Proceeding as in the case of the bimodal distribution, 
one finds that the appropriate generalization of Eq. (|^) 
is given by 



Pn+l(x) 



dl p{l)P n {x — I) 



(36) 



with the same initial condition, Eq. (|4|), as earlier. The 
distribution of the minimal length is then given by 
P™ ln (x) — —dP n (x)/dx and the average minimal length 

is (xr n ) = j °°p n (x)dx. 

Defining P n {x) = G^(x), we recast Eq. (56|) into 



G n+1 (x) = / dlp(l)Gi(x-l). 
Jo 



(37) 



For any given n, G n (x) — > 1 as x — * — oo and G n {x) — > 
as x — * oo. Substituting G n (x) = 1 — g n (x) in Eq. ( |37| ) 
and using normalization of p(l) one finds, 



9n+l(x) 



dyp(x-y)[2g n (y)-gl(y)], (38) 



where g n (x) — ► as x — » — oo and g n (x) — + 1 as x — * +oo. 
We analyze the above equation in the 'backward' tail re- 
gion, i.e., when x — > — oo. In this limit, one can neglect 
the nonlinear term in g n inside the integral in Eq. d3g|). 
The resulting linear equation admits a traveling wave so- 
lution of the form, g n {x) = e x ^ x ~ Vxn \ with 



vx = -\\n 



dz 



,-Ai 



(39) 



For generic length distributions p(l), this function has 
a maximum at A = A* and by the general velocity se- 
lection principle, this maximum velocity is selected, i.e., 
Vmin — v x » ■ 

The question we are interested in is for which class of 
distributions p(l), an unbinding transition can occur for 
the minimal path. Such a transition will occur if the ve- 
locity v\ in Eq. ( |39|) ceases to have a unique maximum. 
By analyzing Eq. (|39[) one sees that this can happen only 
if p{t) has a nonzero delta- function weight at / = 0, i.e., 
when p(l) = a5(l) + (1 - a)f(l) with < a < 1 and /(/) 
is some positive function with dlf(l) = 1. The un- 
binding transition occurs as the parameter a is varied. 
Note that the positivity condition of velocity in Eq. (|39|) 
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demands that 2a < 1. Thus the critical point always 
occurs at a — 1/2. The bimodal distribution considered 
in the previous sections is a special case of this class of 
distributions with a = 1 — p and f(l) = 5(1 — 1). 

For a < 1/2, the average minimal length increases lin- 
early with n for large n with the velocity v\* obtained 
from Eq. (|39|). For a > 1/2, the function P n (x) reaches a 
stationary distribution for large n and hence the (£™ m ) 
saturates to a nonzero constant as n — > oo. What hap- 
pens at the critical point a = 1/2 for generic distributions 
/(O's? For a = 1/2, the Eq. @ reduces to, 



1 



1 



Gn+iW = o G »W + o dlf(l)G 2 n (x - 0. (40) 



Detailed analysis of the above equation reveals a rather 
remarkable universal property at the critical point. It 
turns out that for generic f(l), there are only two pos- 
sible behaviors of (£™ m ) = / °° dxG^(x), depending on 
whether the function /(/) has a gap or not near / = 0. If 
f(l) does not have a gap at I = 0, it turns out that as 
ri — * oo, G n (x) in Eq. ( ^p[ ) tends to a stationary distri- 
bution and hence (£™ m ) — > const. In the opposite case, 
when f(l) has a gap A near I = 0, it turns out that as 
n — > oo, 



A 
ln2 



In In n. 



(41) 



The bimodal distribution corresponds to the special case 
A = 1. 

These general results are best demonstrated by spe- 
cific examples which can be worked out explicitly. Let us 
first consider an example which is gapless, such as the ex- 
ponential distribution, f(l) = e~ l . Remarkably, Eq. ( f4(i| ) 
can be solved exactly in this case. It is easy to verify that 
Eq. ( [i"o| ) admits an exponential solution, G n (x) — A n e~ x , 
where A n+1 = (1 + A 2 n )/2 with A = 0. Asrn oo, A n 
approaches to 1 (more precisely, A n ~ 1 — 2/n). There- 
fore, (L™ in ) = / dxG 2 n {x) = A 2 J2 approaches to 1/2 as 
n — > oo. Next we consider a distribution with gap such 
as f(l) — jS(l — 1) + jS(l — 2). This case also can be 
worked out explicitly following the same steps as used 
for the bimodal case. It turns out that essentially all the 
steps are identical to the bimodal case, except that n in 
Eq. ( p^ ) gets replaced by n/2. But this does not change 
the leading asymptotic behavior for large n which is still 
given by Eq. (p6|). 

We now turn to the maximal path(s) for an arbitrary 
distribution p(l). As in the case of the bimodal distribu- 
tion, R n (x) = Prob(i™ ax < x) satisfies the same equa- 
tion (pi) as P n (x), i.e., 



Rn+l(x) 



dlp(l)R n {x - I) 



(42) 



the only difference is that the initial condition, Eq. 
should be replaced by Eq. (p8|). Substituting R n (x) 
[1 - s n (x)] 2 , we get 



(x) = / dy p(x - y) [2s n (y) - s 2 n (y)] , (43) 



with the boundary conditions s n (x) — > 1 as x — > — oo 
and s n (x) — > as x — > oo. We now have to analyze the 
'forward' tail of the front, i.e., the behavior in the x — > oo 
limit. We first consider bounded distributions p(/)'s with 
an upper cut-off A. In that case, the lower limit of the 
integration in Eq. ( |43| ) becomes x — A. In the tail region 
where x 3> A, one can again neglect the nonlinear term 
in s n inside the integral in Eq. (^3|) and the resulting lin- 
ear equation admits a traveling wave solution of the form 



n{x Hpn) where 



-log 
A* 



2 / dze' xz p(z) 
o 



(44) 



For generic distribution, has a unique minimum at 
p = p* and via the front selection principle, this mini- 
mum velocity is chosen and one gets, i> ma x = ty* . Com- 
paring Eqs. (^9|) and (Q) it becomes evident that for 
the bounded symmetric distributions, i.e., when p(l) = 
p(A — I), one gets the general relation, 



A. 



(45) 



For unbounded distributions, however, one has to be 
careful since it is not obvious that one can neglect the 
nonlinear term in s n inside the integral in Eq. (|4^). 
Equation (Q) for the maximal velocity still remains valid 
as long as p(l) decays with I exponentially or faster. To 
see this explicitly in an example, let us consider the ex- 
ponential distribution, p(l) = e~ l . In this case, one can 



transform the the integral equation (43) into the follow- 
ing difference-differential equation, 



ds n+ i(x) 
dx 



2s n (x) - s 2 n (x) 



(x). 



(46) 



In this differential form, it is clear that for large x, one 
can neglect the nonlinear term on the right-hand side of 
Eq. ([46]). By inserting s n (x) ~ e ~i 1 \ x ~ v ^ n ) into resulting 
linear equation, one gets 



p 



i 



(47) 



Note that the same formula is obtained if one substitutes 
p{l) = e~ L directly into the general expression ([44|). Mini- 
mizing with respect to p, one gets u max - The correspond- 
ing minimal velocity w m i n for the exponential distribution 
can be obtained from the general formula in Eq. (39) with 
p(z) = e~ z and then maximizing with respect to A. So, 
for the exponential distribution we finally get, 



0.231961. 



and 



2.678346 . 



(48) 



Thus for unbounded distributions, there does not seem to 
exist any simple relation between minimal and maximal 
velocities. 
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VI. CONCLUSION 

In this work, we have shown that the length of the min- 
imal path on a Cayley tree exhibits an unbinding phase 
transition from a localized phase to a moving phase. This 
phase transition is driven by the parameter p of the bi- 
modal distribution. In the localized phase (p < p c ), 
the minimal length distribution approaches a stationary 
depth- independent form P mm (a;). This distribution van- 
ishes extremely sharply, as the double exponential, in the 
large x limit. In the moving phase (p > p c ), the mini- 
mal length distribution P™ ln (x) approaches a traveling 
wave form with the front velocity relaxing algebraically, 
Wmin = v(p) + 3/(2A*n) + O (n,- 3 / 2 ) , in the large n limit. 
Specifically, P^ Lltt (x) approaches the solitary traveling 
wave V(y) -V{y + 1), with y = x - (£™ in ), the front 
position (L™ ln ) given by Eq. (|l4|), and V(y) being a so- 
lution of Eq. (||). In the critical case, p = p c , the mini- 
mal length distribution approaches the solitary traveling 
wave, although the front propagates extremely slowly, as 
an iterated logarithm. Given that the distribution of the 
minimal length is the soliton of a finite width, the varia- 
tion of the minimal length is finite. In other words, £™ n 
should not vary much from sample to sample. 

We have also studied the complementary problem of 
the statistics of the maximal path. We have found that 
the maximal length always grows linearly with the height, 
with the velocity w max again determined via the front se- 
lection mechanism. For the bimodal length distribution, 
the minimal and maximal velocities are connected via 
the duality relation, v m - m (p) + u max (l — p) — 1, which 
admits extensions to finite-dimensional situations and to 
arbitrary bounded symmetric length distributions. 

We have found that the tail of the minimal length dis- 
tribution is a double exponential, in contrast with a sim- 
ple exponential that occurs in most traveling wave prob- 
lems ]l6| |. This is not very surprising since similar tails 
were found in statistics of extremes |22j. Note, however, 
that our exact results exemplify extreme value statistics 
for correlated random variables (see |9| for a recent re- 
view), while classical results [E2| correspond to the case 
of independent identically distributed random variables. 

PLK acknowledges support from NSF, ARO, and CNRS. 
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Pn(x) 




X 

FIG. 1. The propagating front for the cumulative distribution P n {x) of the minimal length for the bimodal distribution with 
p = 0.8. 



Rn(x) 




X 

FIG. 2. The propagating front for the cumulative distribution R n (x) of the maximal path for the bimodal distribution with 
p = 0.2. 
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